RKKY Interaction and Intrinsic Frustration in Non-Fermi Liquid Metals 
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We study RKKY interaction in non-Fermi liquid metals. We find that the RKKY interaction 
mediated by non-Fermi liquid metals can be of much longer range than for a Fermi liquid. The 
oscillatory nature of RKKY interaction thus becomes more important in non-Fermi liquids, and 
gives rise to enhanced frustration when the spins form a lattice. Frustration suppresses the magnetic 
ordering temperature of the lattice spin system. Furthermore, we find that the spin system with 
longer range RKKY interaction can be described by the Brazovskii model, where the ordering 
wavevector lies on a shell with constant radius. Strong fluctuations in such model lead to a first- 
order phase transition and/or glassy phase. This may explain some recent experiments where glassy 
behavior was observed in stoichiometric heavy fermion material close to a ferromagnetic quantum 
critical point. 



Introduction: When magnetic moments are placed 
in a metal, the conduction electrons mediate an indirect 
interaction between these moments. Such a long rang in- 
teraction is called the Ruderman-Kittel-Kasuya-Yosida 
(RKKY) interaction. RKKY interaction plays crucial 
roles in, e.g., heavy fermions, diluted magnetic semicon- 
ductors, graphene. The usual derivation of the RKKY in- 
teraction is based on the assumption that the conduction 
electrons form a Landau Fermi liquid. However many 
strongly correlated electron systems show non-Fermi liq- 
uid behavior, e.g. cuprates, heavy fermions, pnictides. 
The question we ask here is what is the form of RKKY 
interaction in a non- Fermi liquid metal, and what are the 
consequencies. 

Of particular interest are heavy fermion systems, where 
local moements couple to the conduction electrons. The 
Doniach phase diagram with competing Kondo coupling 
and RKKY interaction has been the paradigm for heavy 
fermions for decades [T]. In the last few years, as exper- 
imental results accumulate, there is a growing necessity 
to go beyond the Doniach phase diagram. Frustration 
or the quantum zero point energy has been proposed as 
a new dimension in the global phase diagram of heavy 
fermions [2HS] ■ One obvious origin of frustration is frus- 
tration of lattice structure itself. However such geometric 
frustration is not universally observed in heavy fermion 
materials. Here we propose that the non- Fermi liquid 
nature of conduction electrons in the Kondo liquid phase 
leads to intrinsic frustration for the localized spin de- 
grees of freedom. This provides a more universal source 
of frustration. 

Our approach is based on the idea of quantum critical- 
ity and non-Fermi liquid (NFL) behavior. The standard 
picture is that the critical fluctuations near a quantum 
critical point (QCP) leads to NFL behavior. Here we 
depart from this picture by starting with the assump- 
tion that in a certain range of the parameter space, the 
itinerant electrons form a NFL state. We then proceed 
to study its consequences on other degrees of freedom, 
e.g. the localized spins. Focusing on the regime with 
small Kondo coupling, i.e. a small Fermi surface, we find 
that the magnetic transition temperature will be reduced 
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FIG. 1: Schematic electronic and magnetic phase diagrams. 
Distance dependence of the RKKY interaction is shown in the 
insets. In the non-Fermi liquid region, RKKY interaction is of 
longer range, leading to frustration. The magnetic transition 
temperature decreases with increasing frustration, and new 
phases can emerge near the QCP. 



by the frustration resulting from longer-range RKKY in- 
teraction produced by NFL itinerant electrons. Further- 
more, we find that the putative ferromagnetic (FM) QCP 
may be replaced by a first-order phase transition or a 
glassy phase [718] (see Fig. [I]). 

Formalism: We start with the Kondo lattice model, 
H = H c + H K . Here H c is the conduction elec- 
tron Hamiltonian, and ususally only the hopping term 
is included. The Kondo coupling between conduction 
electrons and localized spins is of the form, Hk — 



Ik 



J^iapS* ' c ia a a0 c i0- We depart from the usual 
approach by considering the conduction electrons to be 
strongly interacting themselves, i.e. He = Hq + H ( ™ ' . 
One way to motivate this is to consider the phenomeno- 
logical two fluid model [5HT2]. In many heavy fermion 
systems, below the coherence temperature T* , the ex- 
perimental results can be understood in terms of the 
two fluid model, with one component the itinerant heavy 
electrons, and the other component local moments. The 
heavy electron Kondo liquid is not a simple Fermi liquid, 
e.g. the specific heat is logarithmically enhanced at low 
temperature. One has a model of interacting itinerant 
electrons coupled with localized spins. 

The itinerant electrons induce a RKKY type interac- 



tion among the localized moments: 



tj \ s jab nfl nb 

-nRKKY — 2_^ J ij b i ^j ■ 



(1) 



Here the coupling J^ = —^xtj [T5HT5] . is determined 



Jk , ,afc 
4 X-ij 

by the static spin susceptibility of the conduction elec- 
trons 



Y ab 



([ S Q (n,i), S 6 (r J ,0)])e-" t dt, (2) 



with the electron spin s a (n) = J2 a p c \ a a a$ c iP and 
rj = + . If the conduction electrons are in the param- 
agnetic state, the spin susceptibility is isotropic and di- 
agonal, i.e. Xij = x( r ij)d ab - F° r Fermi liquids, the spin 
susceptibility behaves as x( r ) ~ (l/?* d ) cos(2/«i?r + 0q) 
at long distances, with d the spatial dimension. This 
leads directly to the standard form of the RKKY inter- 
action. The exponent d results from the sharp jump in 
the momentum distribution n(fe), characteristic of Fermi 
liquids. For non-Fermi liquid metals, the RKKY inter- 
action can have qualitatively different behavior. We still 
assume the existence of a Fermi surface, i.e. singularity 
in n(k), thus the spin susceptibility still has 2kp oscilla- 
tion. The exponent can take a different value. Thus we 
have x( r ) ~ {l/r a )cos(2kpr + 8o). More detailed studies 
of the NFL spin susceptibility will be presented below. 

Consider placing a lattice of spins in the non-Fermi 
liquid metal. We focus on the effect of the RKKY in- 
teraction to the spin system, and will not consider the 
competition between Kondo coupling and RKKY in- 
teraction pp. This can be achieved by assuming the 
spins to be classical, or considering only the part of 
the phase diagram with a small Fermi surface. With 
S(q) = (l/N)J2i Sie lq " ri , one has in momentum space, 
H = j: q F(q)S(q)-S(-q), where 



F (?) = ^ E j ^y q ' 



(3) 



r-i#0 



with Ti defined on the lattice. The ordering wavevec- 
tor in the ground state is determined by minimizing the 
function F(q). 

For the conventional three dimensional RKKY interac- 
tion mediated by a Fermi liquid, this problem has been 
studied in |16j . where different phases have been iden- 
tified as the conduction electron density changes. At 
small kpa, where a is the lattice constant, the ground 
state is ferromagnetic. As kpa increases, antiferromag- 
netic phases with different ordering wavevectors appear. 
In the case kpa — ¥ 0, the above summation can be re- 
placed by an integral, and F(q) ~ — xil)- The ordering 
wavevector is thus determined by maximizing the static 
spin susceptibility. 

Now we consider in more detail what is the form of 
the static spin susceptibility in a NFL metal. When 
vertex corrections can be ignored, the spin susceptibil- 
ity can be calculated from the fermion bubble, with 



Xab(q) ~ J a a G(k + q)<j b G(k). When the momentum 
distribution n(k) has a weaker singularity than a jump 
at kp, e.g. a kink, the Friedel oscillation decays faster 
than that of Fermi liquids (see Appendix). Then one ex- 
pects x( r ) an d J(r) to decay faster than that of Fermi 
liquids. An interesting question is whether it is possible 
to have longer range RKKY interactions, which would 
generate the desired frustration among the spins [2-6J. 
We will present two models of NFL metals that can give 
rise to such behavior. 

Longer range RKKY interaction in 1-d: First, 
as a proof of principle that RKKY interaction in a 
strongly interacting electron system can be of longer 
range than in a free system, let us first consider one di- 
mension. In 1-d, RKKY interaction mediated by free 
electrons is of the form J(r) ~ Si(2fci?r) — |, with 
the sine integral function Si(x). At large distance one 
has J(r) ~ cos(2kpr)/r. In momentum space, one has 

x{q) ~ (1/g) l n l 9 ^2fc F I' w ^ tn a max i mum at 1 — 2kp. 

The low energy dynamics of interacting electrons in 1-d 
is described by the Luttinger liquid theory. Due to spin- 
charge separation, the conduction electron Hamiltonain 
can be written as a summation of the two channels |17j , 



H c= E y/<M3«n 2 Q 



g^{d x e a f 



(4) 



with v c and v s the velocity of charge and spin density 
wave respectively. The charge interaction constant g c — 1 
for noninteracting fermions, g c < 1 for repulsive interac- 
tion, and g c > 1 for attractive interaction. We are in- 
terested in the case with repulsive interaction. The spin 
interaction constant g s = 1 in the presence of SU(2) spin 
symmetry The oscillating part of the spin correlation 
function is [T7] 



(s(x,t) ■ s(0,0)) 



cos(2kpx) 



\t + ix/v c \3" \t + ix/v s \ 9s 



(5) 



The RKKY interaction, determined from the static spin 
susceptibility, is of the form 



J(x) 



dr 



cos(2kpx) 



\t + ix/v c \ 9c \T + ix/v a \ gs 



cos(2kpx) 



, (6) 
For g c < 1, g s = 1, the exponent a = g c + g s — 1 < d = 1. 

The RKKY interaction mediated by a Luttinger liquid 
is thus of longer range than that mediated by a non- 
interacting Fermi gas^H 

Spin susceptibility in 2-d: Now we consider two 
dimenional metals. For free electrons, the static spin 
susceptibility reads 



X(q) 



Xo for q < 2k F 

Xo [l - y/l - {2k F /q)A for q > 2k F , 



(7) 



with xo — l/ 71 ") which has a one-sided square-root sin- 
gularity. The RKKY interaction is thus of the form 
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FIG. 2: (a) The static spin susceptibility x(l) as function of 
momentum for Fermi liquid (dashed, black) and the gauge- 
fermion model with a < 1/3 (dotted, blue) and a > 1/3 
(solid, red), (b) F{q) as function of momentum for an- 
gles = 0, 7r/6, 7r/4. The curves for different angles are al- 
most identical. Here the spins form a square lattice, and 
a = 1/2, k F a = 0.2. 



J(r) ~ sin(2k F r)/r 2 W^ For 2-d Fermi liquid, including 
higher order diagrams, there is also a square-root singu- 
larity for q < 2k F , with x(q) = x(2fc F ) + X sins («) [IB], 
and 



X sing («7) - iv/1 - (<?/2M 5 



(8) 



The new singularity gives contribution 5x{i") ~ 
[(2fci?r) cos(2/cpr) — sin(2/c^r)]/r 3 , and will not change 
the long distance behavior of the RKKY interaction. 

A prototype of non-Fermi liquid metal in higher di- 
mensions is the system of 2-d degenerate fermions inter- 
acting via a singular gauge interaction [ITJH2"rj] , where the 
presence of the gauge interaction leads to singular 2k F 
response |23j . The fermion 2k F vertex T<ik F has a power 
law dependence on frequency, with T 2 k F ~ ( — 



Iki 

l 

2N 



1 2k f 



In- 3 N 



The exponent is of the form a — -^ -r 1;r2jv - 2 
O (^) for large N, and a = ^| + O (1) for small N. 
Here the spin index is generalized to take values from 1 
to N. Taking N — 2, one obtains a = 0.35 from the 
large-N expansion, and a = 0.56 in the small-N limit. 

The spin susceptibility is calculated from the polariza- 
tion bubble with vortex corrections, x(g,u>) — n(q,w) = 
jdpdeG{p+q/2,e + Lo/2)G{p-q/2,e-u J /2)[T tp {q,u J )] 2 . 
For a < 1/3, the static spin susceptibility is of the form 



X (q)~X0-C\q-2k F \ 1 - 3 °, 
and for a > 1/3 one has [2"3] 



\q-2k F \^- 1 ' 



(9) 



(10) 



with a singularity at g = 2fej? (see Fig. [2^ a)). Fourier 
transforming to real space, we find 



X(r) 



1 cos(2k F r - O ) 

iJo{qr)qdq~ 



\q-2k F \^ 



r 5/2-3o- 



The exponent a = 5/2 — 3a can be much smaller than 
the space dimension d = 2. 

More generally, for non-Fermi liquid metals, one can 
employ a scaling theory for the susceptibility (see e.g. 



[27128) ). Assuming the existence of a Fermi surface, the 
static spin susceptibility generally has a power law behav- 
ior near q — 2k F , with x(q) ~ \q ~ 2k F \ v . For v < 1/2, 
one has a stronger singularity than the Fermi liquid case, 
and the RKKY interaction is of longer range. 

Longer range RKKY interaction in 2-d: Let us 
now consider the ground state of the spins embedded in 
the 2-d metals with small k F a. For the Fermi liquid case 
(Eqs.( |7|8| ), x(q) increases monotonically with decreasing 
q (see Fig. [2]). The ground state is ferromagnetic. For 
non-Fermi liquid (Eqs.([9]),( 10)), the maximum of x(l) ^ s 
at q = 2k F , indicating an instability of the ferromagnetic 
state. More precisely, one can calculate the interaction 
F(q) by first Fourier transforming x(q) to real space to 
get x( r ): an d then performing the lattice summation in 
Eq. Q. The result for a = 1/2 is shown in Fig. ^b). 
One can see that F(q) has a minimum at q = 2k F . The 
singularity in xil) i s smeared out by the lattice effect. 

Another observation is that F(q) has a very weak de- 
pendence on the direction of momentum. In Fig. [2Fb) , 
F(q) for the three different angles are almost indistin- 
guishable. With the minimum of F(q) at qo — 2k F , the 
ordering wavevector of the lattice spin system lies on a 
shell of radius 2k F . Expanding F(q) around qo, one ob- 
tains the Brazovskii model 1291. 



H = J2[bo + D (\q\ - go) 2 ] S(q) ■ S(-q) 



(12) 



Brazovskii found that the large phase space available for 
fluctuations around a shell of minima leads to a first- 
order phase transition )29) . It has been found experi- 
mentally that putative FM-QCPs are replaced by first 
order transitions at low temperatures in several transi- 
tion metal compounds, e.g. MnSi, ZrZn.2, and heavy 
fermion systems, e.g. UGe 2 , UCoAl, UCoGe (see [30] 
and references therein). It was realized earlier that com- 
peting orders [3_T] as well as fluctuations [32-55] can lead 
to first order quantum phase transitions. Here we find 
a new mechanism where the frustration resulting from 
NFL behavior generates first order transitions. 

A further observation is that the extensive configura- 
tional entropy in the Brazovskii model should lead to 
slow dynamics and glassiness [36H39] P^ Glassy correla- 
tions emerge when the correlation length £ = [D/b) 1 / 2 
becomes of order the modulation length l = 2n/q [37] . 
The parameter b needs to be determined self-consistently. 
Within the large-N approximation, and including a small 
quartic term with coupling u, we have 



b = bn + uT 



(27T) 



Mq), 



(13) 



with the Green's function Q(q) = l/[b + D(q — q ) }. 
The condition £,/Iq ~ 1 then determines the temperature 

where glassy behavior sets in to be T g ~ ^r~ ci-'iog/qoa) > 
with the coefficient c\ of order unity and momentum cut- 
off A ~ aT 1 . We notice that here T g depends logarithm!- 



cally on cutoff instead of the 1/A dependence for the 3-d 
model considered in |37j . 

Glassy spin dynamics was recently observed in the 
heavy fermion system CeFePO [3D]. CeFePO is a layered 
Kondo lattice system, in close proximity to a FM QCP. 
Spin-glass-like freezing was detected in the ac suscepti- 
bility, specific heat and muon-spin relaxation [40 . The 
glass behavior in such a stoichiometric system points to 
new mechanisms that do not reply on external random- 
ness. Our model provides such a possibility (see [4TU46] 
and references therein for earlier attempts to obtain glass 
behavior from frustrated deterministic models). 

Away from QCP: Having identified a glass transi- 
tion near the QCP, we proceed to study the behavior 
of the lattice spin system away from QCP using a ran- 
dom exchange model that is widely used to describe spin 
glasses. Due to the cosin function, the RKKY interac- 
tion changes sign and magnitude with distance. It can 
be well approximated by a random interaction 

T 2 
J V 4 , - 

distribution P{e vj ) = (1/tt)(1 - e^-) 

When the itinerant electrons are away from the QCP, 
there is a crossover to the Fermi liquid behavior at low 
energy, or equivalently long distance, where the RKKY 
interaction is substantially reduced. We will assume for 
simplicity that the RKKY interaction can be neglected 
beyond a crossover scale rpL- Then the exchange inter- 
action is of the form 



where e^ is a random variable with cosine 

-1/2 



■~>ij — 







for 
for 



> r FL . 



(14) 



When rpL becomes of order the lattice constant, only 
the nearest neighbor interactions survive, i.e. H = 
J^2<u> Sj • Sj, and the spins are magnetically ordered. 
As rpL increases, the ordering temperature will be re- 
duced by frustration. 

A simpler model that illustrates essentially the same 
effect of suppression of ordering by frustration is the 
Sherrington-Kirkpatrick model [51 52]- Consider here 
ferromagnetic ordering. We start with a mean field type 



Hamiltonian H 



-JoE 



m 



Sj, with Jo > 0, and 



each spin interacts with z neighouring spins. The spins 
order ferromagnetically below the transition temperature 
T c (0) = J S(S + l)/6, with J = zJ . This correpondcs 
to the case far away from the QCP. 

Then we add to the above mean field ferromagnetic 
model random exchange interactions to model the frus- 
tration effect when approaching a QCP. The new Hamil- 
tonian can be written as H = — Y^(a) JijSi ' &ji where 
the interaction J^- is distributed according to P(Jij) = 
{ \ji o)2 \ |5T153] . This model is readily 

solved by the replica technique [52153) . and the transi- 
tion temperature to ferromagnetism is reduced by the 
random interactions, with the result |52|54) 



V27rJ ; 



exp 



T = T (o) 




(15) 
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FIG. 3: Ferromagmetic transition temperature as function of 
range of random exchange interaction. 



where we have defined J = z 1 / 2 J . 

We fix the mean field ordering temperature in the ab- 
sence of random exchange interaction T c and the vari- 
ance of the random distribution J, so that z is a measure 
of t he r ange of random exchange interaction, i.e. tfl in 
Eq.Q. We can define z c = {S{S + l)/3) J^/J 2 , and 
write T r in the form 



T c = Tj°) 



1- 



(16) 



which is plotted in Fig. [3j One can see that with in- 
creasing range of random exchange interaction, the FM 
ordering temperature decreases. This then translates to 
the picture that when approaching the QCP, as RKKY 
interaction becomes of longer range, magnetic ordering 
is suppressed (see Fig. [I]). 

Conclusions: We have studied the RKKY interac- 
tion in non-Fermi liquid metals. The basic picture we find 
is summarized in Fig. [I] In the non-Fermi liquid phase, 
when including vertex corrections, the RKKY interaction 
can be of longer range than in a Fermi liquid. Longer 
range RKKY interaction leads to frustration for the lat- 
tice spin system placed in such a NFL metal. Magnetic 
ordering will be suppressed by frustration, and novel be- 
havior may emerge near the putative QCP. In particu- 
lar, the continuous second-order phase transitions may 
be replaced by first-order transtions. Glassy dynamics 
may occur near the QCP without invoking disorder. One 
candidate material for such glass behavior is the heavy 
fermion system CeFePO. We focused here on FM QCP. 
One can also generalize the whole procedure to AFM 
QCP by increasing kpa. Another interesting question 
is the competition between the Kondo coupling and the 
longer range RKKY interaction. 
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Appendix: Priedel oscillation 



When n(k) has a jump at Uf, one has 



We consider here the Friedel oscillation for a non- Fermi 
liquid metal. When the momentum distribution n(fe) has 
a weaker singularity than a jump at Uf, e.g. a kink, the 
Friedel oscillation decays faster than that of Fermi liq- 
uids. Consider for simplicity d — 2, and Fourier trans- 
form gives 



n(r) 



fcJi(fcr) 



n(k) 



k F +S 



k F -6 



cos(kpr — 37r/4) 
{k F rf/ 2 



(18) 



When n(k) has a higher order singularity, further partial 
integration leads to 



n(r) ~ / dkkn(k)Jo(kr) , 



(17) 



with Jo the Bessel function. Contribution from the sin- 
gularity of n(k) can be obtained by partial integration. 



n(r) 



kJ%(kr) dn(k) 
r 2 dk 



k F +8 



k F — 8 



cos(k F r - 5tt/4) 



which decays faster than the Fermi liquid result. 
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